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^! , Abstract. In the present paper, we introduce a two-component model of the Galactic 
disk to investigate its chemical evolution. The formation of the thick and thin disks 
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1. Introduction 

Since the existence of the thick disk of our Galaxy was confirmed by Gilmore & Reid 
(1983) more than ten years ago, it has been generally accepted that a complete de- 
scription of the thick disk, such as its scale length, scale height, density normalization, 
metallicity and kinematical properties, is a necessary step towards understanding the 
Galaxy formation, halo collapse, disk dynamical and chemical evolution. Unfortunately, 
the characteristics of this population remain controversial, especially for its density pro- 
file. Several attempts have been made to deduce this parameter by remote star counts 
and field star survey ( Buser & Kaeser 1985, Gilmore, Reid & Hewett 1985, Reid & 
Majewski 1993, Buser, Rong & Karaali 1998). But the results are still quite uncertain. 
This is partly due to the lack of complete sample of thick disk stars since its members 
cannot be easily recognized from that of the thin disk and/or the halo in most observable 
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distributions. Moreover, the determination of thick disk characteristics requires large star 
samples in various directions well distributed in both the longitude and latitude (Robin 
et al. 1996), which cannot be obtained easily at the present time. 

Recently, the study of chemical evolution of the Galactic disk has been proven to be 
a powerful tool to explore the formation and evolution of our Galaxy. Numerous models 
have been detailedly put forward (Matteucci & Francois 1989, Ferrini et al. 1994, Gio- 
vagnoli & Tosi 1995, Prantzos & Aubert 1995, Timmcs et al. 1995, Carig 1996, Pilyugin 
& Edmunds 1996a, Chiappini, Matteucci & Gratton 1997, Allen et al. 1998, Thon & 
Meusinger 1998, Prantzos & Silk 1998). Among them, Chiappini, Matteucci & Gratton 
(1997 hereafter CMG97) were the first to take into account the effect of thick disk. They 
assumed that there are two main accretion episodes. The first is responsible for the for- 
mation of the thick disk, and the second, delayed relative to the first, forms the thin 
disk. The predictions of their best-fitting model are in good agreement not only with the 
observed metallicity distribution , but with the observed number of very low metallic- 
ity stars (Rocha-Pinto & Maciel 1996, hereafter RM96). This enlightens us to do more 
detailed analyze of the disk evolution based on the new chemical constraints. 

In the present paper, the two-component model for the Galactic disk evolution (such 
as CMG97) is adopted, in which the local surface density of the thick disk at the present 
time is chose to be one of the free parameters. The infall rate is assumed to be in a 
Gaussian form instead of an exponentially decreasing one. The quantitative compari- 
son between model predictions and the observations, i.e., the new G-dwarf metallicity 
distribution obtained by Hou et al (1998), is used for the x 2 -test of the best-fitting model. 

The outlines are as follows. In section 2, we present brief description of observa- 
tional constraints up to now, of which the most important is the G-dwarf metallicity 
distribution. Section 3 is the model and its main ingredients. In section 4, we present 
best-fittings of four different models, which are closed-box, one-component, pre-thin and 
post-thin models respectively, to the observations. Discussions of the models are also 
included in section 4. Our conclusions are shown in the last section. 

2. Observational constraints 

A successful model of the chemical evolution of the galactic disk should reproduce the 
main observational features of both the solar neighborhood and whole disk. Our set of 
constraints includes: 

(1) G-dwarf metallicity distribution in the solar neighbourhood (Hou et al 1998); 

(2) radial abundance gradients at present time; 

(3) age- metallicity relation (AMR); 
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(4) the correlation between [O/Fe] and [Fe/H]; 

(5) radial profiles for the gas surface density; and 

(6) the variations of Star Formation Rate (SFR) across the disk. 

The first one is selected as the observational constraint to quantitatively estimate 
the best-fitting model in this paper, since G dwarfs cover the whole life of the Galactic 
disk and its metallicity distribution can reflect the local chemical enrichment history 
The others are used for the comparisons between the best-fit model predictions and 
observations. 

2.1. G- dwarf metallicity distribution 

The metallicity distribution of G dwarfs in the solar neighborhood is one of the most 
important constraints on the chemical evolution of the Galactic disk. Since G dwarfs 
have lifetimes comparable to the estimated age of the Galaxy, they represent a sample 
which has never been depleted by stellar evolution, accumulating since the first episodes 
of low-mass star formation. Therefore, a complete sample of these stars in the solar 
neighborhood carries memory of the local star formation and chemical enrichment history 

Pagel & Patchett (1975) derived a cumulative G-dwarf metallicity distribution, based 
on a volume-limited sample of 132 G dwarfs within about 25 pc of the Sun. Pagel (1989) 
revised previous data of Pagel & Patchett (1975) by means of a new calibration between 
the ultraviolet excess 5 (U-B) and [Fe/H]. Later, Rana (1991) and Sommer-Larsen (1991) 
independently revised the distribution of Pagel (1989), taking into account the dynamical 
heating effect on the observed distribution. RM96 derived a G-dwarf metallicity distribu- 
tion in the solar neighborhood, using uvby photometry and up-to-date parallaxes. RM96 
introduced a chemical criterion, according to which all sters have [Fe/H] < -1.2 are con- 
sidered to be halo members and excluded from the final sample. The distribution of 
RM96 comprises 287 G dwarfs within 25 pc from the Sun and differs from the classic one 
by having a prominent single peak around [Fe/H] = - 0.2 (sec RM96 for details). 

Recently, Hou et al (1998) collected a new, enlarged sample of G dwarfs within 25pc 
from the Sun. The stars are selected from the third Catalogue of Nearby Stars (Gliese & 
Jahrciss 1991). The uvby data are taken from the catalogues of Olsen (1993) and Hauck 
& Mcrmilliod (1990). No chemical criterion was introduced in Hou et al (1998) since 
observational evidences showed that, in the metallicity interval -1.5 < [Fc/H] < -1.0, the 
fraction of thick disk stars in the solar neighbourhood appears to be as high as 60% 
(Nissen & Schuster 1997). This is one of the main differences between the distribution 
of Hou et al (1998) and that of RM96. The final sample contains 382 G dwarfs with 
photometric data. The adopted metallicity calibration and kinetic correction is the same 
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Table 1. The obtained G-dwarf metallicity distribution in the solar neighbourhood taken 
from Hou et al (1998) 
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1.00 
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as that of RM96. Following Pagel (1989) and RM96, Hou et al (1998) have also corrected 
the distribution for observational errors and cosmic scatter. The results of Hou et al (1998) 
are shown in Table 1, in which the sample is divided into 19 bins. The first column of 
Table 1 is the metallicity range in each bin. The raw distribution (A7Vo)j are presented 
in column 2. The third is the weight factor for the scale height correction according 
to Sommer-Larson (1991). The fourth presences 5(ANo)i, the correction factors for the 
observational errors and cosmic scatter. The obtained relative distribution is given in the 
last column, where N = X^it ^J^ + 3{AN a )i] is the total number of G dwarfs after 
correction. In the last column, the second to fourth bins are grouped to yield a mean value 
for the distribution, which is similar with the method used in RM96. This procedure is 
also used for the fifth to seventh bins. Table 1 shows that the resulted distribution of Hou 
et al (1998) differs from that of RM96 by having a larger width and smaller amplitude of 
the single peak. Moreover, the metal-poor tail of the new distribution (Hou et al 1998) 
extends to [Fe/H] =-1.5. 
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2.2. Abundance gradients 

Furthermore, radial metallicity variations of the interstellar medium (ISM) can constrain 
models of Galaxy formation and chemical evolution. From extensive studies of optical 
emission lines in HII regions, Shaver et al. (1983) derived an oxygen abundance gradient 
of the order of -0.07 dex/kpe. Afflcrbach et al.(1996) have indirectly deduced a similar 
result by measuring electron temperature variations in a set of ultra-compact HII regions. 
A relatively flatter gradient has been obtained by Vilchez & Esteban(1996) for the outer 
Galaxy, using spectroscopic observations of a sample of HII regions towards the Galactic 
anti-center. 

The recent radial profile of oxygen in the Galaxy can also be traced by observations 
of B-type stars, with main sequence ages less then 1.0 Gyr. A series of medium- to high- 
resolution spectroscopic observations of early B-type main-sequence objects have been 
published (Smartt et al. 1997 and references there in). Using this homogeneous sam- 
ple, Smartt et al.(1997) derived an oxygen abundance gradient of -0.07 ± 0.01 dex/kpc 
between Galactocentric distance of 6 kpc < r < 18 kpc, which is in good agreement 
with nebular studies. Gummersbach et al (1998) determines the stellar parameters and 
abundances of several element for 16 early B main-sequence stars in Galactocentric dis- 
tance 5.0 kpc < r < 14.0 kpc by reanalyzing and extending the observations of Kaufer 
et al (1994). An oxygen abundance gradient -0.07 ± 0.01 dex/kpc is derived, typical for 
normal spiral galaxies of similar Hubble type. 

2.3. Others 

Just as mentioned above, other chemical constraints should be taken into account at 
the same time, such as the age-metallicity relation, the correlation between [O/Fc] and 
[Fe/H] for field stars as well as the radial profiles of gas surface density and SFR at the 
present time for the whole disk. 

Twarog (1980) obtained the first AMR for the local disk stars, with the stellar ages 
from the theoretical isochrones. The same sample has been reanalyzed by Carlberg et 
al. (1985) with different results, due to the revision of the isochrones as well as the 
calibration of abundances. A more accurate AMR was obtained by Edvardsson et al. 
(1993). They have derived abundances of 13 different elements, such as O, Fe, Si, Ba 
etc., as well as individual photometric ages, for 189 nearby field F and G dwarfs. Their 
abundance analysis was made with theoretical LTE model atmospheres, based on the 
extensive high resolution, high S/N, spectroscopic observations of carefully selected field 
stars. The resulted AMR of Edvardsson et al. (1993) was used at the present study. 
However, this AMR does not constitute a tight constraints of the chemical model, since 
there is a considerable scatter. Moreover, the results of the survey of Edvardsson et 
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al. (1993) concerning O vs. Fe relationships for field stars are used in this study. As 
for metal-poor stars, the correlation between [O/Fe] and [Fe/H] are taken from Barbuy 
(1988). 

The radial Galactic profiles of atomic and molecular hydrogen are discussed in Lacey 
& Fall (1985). An updated discussion is given in Dame (1993). Inside the solar circle, 
the molecular and atomic gas are found in roughly equal amounts. However, the surface 
density of atomic hydrogen, which seems to be constant from Akpc to 15fcpc, dominates 
the gas profiles outside the solar circle. The radial distribution of the sum of atomic and 
molecular hydrogen given in Dame (1993) is adopted in this paper. 

The radial distribution of the present SFR in the Galaxy are taken from Gusten & 
Merger (1983), Lync et al (1983) and Guibert et al (1978). Data are based on several 
tracers of star formation: Lyman continuum photons from HII regions, pulsars and su- 
pernova remnants. It is normalized to the present SFR in the solar neighbourhood (as in 
Lacey & Fall 1985) since the absolute values depend on poorly known conversion factors. 

3. The model 

It is assumed that the Galactic disk is sheet-like, which originates and grows only from 
the infall of primordial gas. The disk is considered as a system of independent rings with 
lkpc wide for each. No radial inflows or outflows are considered and the center of each ring 
locates at its median Galactocentric radius. The ring centered at Galactocentric distance 
r Q = 8.5 kpc is labeled as the solar neighbourhood. The age of the disk is adopted to be 
13.0Gyr (Rana 1991). 

3.1. Basic equations and nucleosynthesis 

The instantaneous-recycling approximation (IRA) is relaxed, but instantaneous mixing 
of the gas with the stellar ejecta is assumed, i.e., the gas is characterized by a unique 
composition at each epoch of time. We solve numerically the classical set of equations of 
Galactic chemical evolution (Tinsley 1980, Pagel 1997) as 
dT, tot (r,t) 



dt 



= f(r,t), (1) 



' ' =-i>(r,t)+ (m-m r )^(r,t-r m )^(m)dm + /(r,t), (2) 



d[Zi(r,t)Z gas (r,t)] 



')S (r t)\ f 

^ dt 9asK 1 n = -Zi(r, i)V(r, t) + J m yi , m ^(r, t - T m )4>(m) dm + Z ilf f(r, tp) 

where T, tot (r, t) and S gas (r, t) are the total and gas surface density respectively in the ring 
centered at Galactocentric distance r at evolution time t; f(r,t) is often called the infall 
or accretion rate; ")p(r,t) is the star formation rate (SFR) and 4>{m) is the initial mass 
function (IMF); m r and T m are the remnant mass and the lifetime of a star of initial 
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mass m, respectively, and m t is the corresponding initial mass of a star whose main- 
sequence lifetime r m equates to evolution time t (the turnoff mass). Here the mass range 
of IMF is taken from O.1M to 100M Q . The mass of element i in the gas evolves via star 
formation (putting metals from the ISM into stars), ejection, and gas inflows, according 
to equation (3), where t/j ;TO is the stellar yield of element i, i.e., the mass fraction of a 
star of initial mass to that is converted to element i and ejected, and Zij is the mass 
abundance of element i in the infalling gas, which is assumed to be piromordial in this 
study: Zoj — Zp e j — 0. It should be emphasized that the second terms in the right 
hand of equations (2) and (3) also include the contribution of Type la supernovas (Type 
la SNs), which is detailcdly presented in Matteucci & Greggio (1986). The constant A in 
equation (9) of Matteucci & Greggio (1986) describes the fraction of systems with total 
mass in appropriate range, which eventually succeed in giving rise to a Type la SN event, 
and in this study, it is fixed by requiring to present best-fit to the metal-rich tail of the 
G-dwarf metallicity distribution in the solar neighbourhood. 

It is also assumed that every star ejects its envelope just after leaving the main 
sequence. The adopted relation between main-sequence lifetimes r m ( in units of Gyr) 
and stellar initial mass to (in units of M Q ) is ( Larson 1974): 

logm = 1.983 - 1.054y/(logT m + 2.52). (4) 

For the sake of simplicity, we assume that, except for Type la SNs, any star evolves as 
a single star even if it is the member of a binary system. All massive stars (to > 9Mq) 
explode as type II supernovas (Type II SNs), leaving behind a neutron star of mass 
m,R = O.5M (Prantzos & Silk 1998). The final stage of the intermediate /low mass 
stars (M < 9M Q ) is white dwarfs, and the final-initial mass relation is taken from 
Weidcmann(1984). Type la SNs are thought to originate from carbon deflagration in 
C-0 white dwarfs in binary systems. The method included the contribution of Type la 
SNs is the same as Matteucci & Greggio (1986). 

In this paper, we only consider the evolution of iron and oxygen. The oxygen and iron 
production for Type II SNs and Type la SNs are taken from Woosley & Weaver (1995) 
and Woosley (1997), respectively. Recently, using the evolutionary tracks of Geneva group 
up to the early asymptotic giant branch (AGB) in combination with a synthetic thermal- 
pulsing AGB model, van den Hock & Groenewegen (1997) calculated in detail the chem- 
ical evolution and yields of six elements up to the end of AGB. Their results showed that 
the low-mass stars (to < 3M Q ) produce small amounts of oxygen, yet the intermediate 
mass stars (3M < to < 8M ) destroy the initial oxygen through Hot Bottom Burning 
(HBB). Therefore, it is reasonable in this paper to neglect the oxygen production by 
intermediate/low mass stars compared with that of massive stars. 
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3.2. The infall rate 

Currently popular models of the galaxy formation are semi-analytic models within the 
framework of the hierarchical structure formation paradigm (White & Rees 1978, White 
& Frenck 1991, Wechsler et al 1998, Mo et al 1998), which allow one to model the 
astrophysical processes involved in galaxy formation in a simplified but physical way 
(Kauffmann et al 1993, Somerville & Primack 1998, Kauffmann et al 1998, Primack et al 
1998). These models are in good agreement with a broad range of local galaxy observa- 
tions, including the correlation between luminosity and circular velocity for spirals ( the 
Tully-Fisher relation), the B-band luminosity function, cold gas contents, metallicities, 
and colors. These models postulate that the formation of galaxy is mainly regulated by 
gas cooling, dissipation, star formation, and supernova feedback. However, for the pur- 
pose of simplicity, Galactic chemical evolution models assume that the integral effect of 
these processes can be represented by that of an infall rate, which is a function of evolu- 
tion time and Galactocentric distance. The form of the infall rate is adjusted to satisfy 
the constraint of G-dwarf metallicity distribution in the solar neighbourhood. Although 
an exponentially decreasing infall rate is widely used, we adopt here a Gaussian form 
for the infall rate. The physical motivation for such a choice is that because of its small 
initial surface density, the local disk initially accretes a small amount of the surrounding 
gas; as the disk mass and gravitational potential build up, the accretion rate gradually 
increases, but starts decreasing when the gas reservoir is depleted (Prantzos & Silk 1998). 

Numerous models for the formation of thick disk have been put forward since the 
confirmation of its existence by Gilmore & Reid (1983) (see Majewski 1993 for details). 
The models fall into either " top-down" scenarios (the pre-thin model) , where the forma- 
tion of the thick disk precedes that of the thin disk, or the " bottom- up" scenarios (the 
post-thin model) , where the thick disk is the result of some action on or by the thin disk. 
In the pre-thin disk model, the formation of the thick disk is a transitional phase during 
the general contraction of the Galaxy. This model views the thick disk as a dissipative, 
rotational-supported structure, and the halo as non-dissipative and supported by the 
kinetic pressure provided by large, anisotropic velocity dispersions. The post-thin model 
resorts to formation of the thick disk after the gas has completely collapsed into a thin 
disk. Possible physical processes are: (1) Secular kinematic diffusion of thin disk stars 
(Norris 1987); (2) Violent thin disk heating by the accretion of a satellite galaxy (Quinn 
et al 1993). The required events must not occur too late in the disk life time so that the 
gas can cool again and form stars in the thin disk; (3) Halo response to disk potential 
(Gilmore & Reid 1983). The post-thin model leaves two important observational signa- 
tures. First, the thick disk is a separate population distinct from the thin disk and the 
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halo. Second, no gradient can be generated in the thick disk by the events, although a 
pre-existing gradient may survive the merger. 

In this study, both the pre-thin and post-thin model for the formation of Galactic 
disk are considered. Following CMG97, we assume that there are two main infall episodes 
in both cases. The rate of mass accretion (in unit of M^pc^Gyr -1 ) in each ring could 
be expressed as 

/M) = A^L e -(t-r t f^ + J^L e -(t-T d r/2*i (5) 

V27TCr t \/2lTOd 

where Tt and (in units of Gyr) are the maximum infall time of the thick and thin disk 
respectively, at and <Jd (in units of Gyr) are the corresponding half-width. It is assumed 
that at ~ r t and ad ~ Td, i.e., the value of the half- width can just varies a little around 
that of the formation time-scale (Prantzos & Silk 1998). 

For the pre-thin model, the first infall episode forms the thick disk, which originated 
from a fast dissipative collapse such as that suggested by Eggen, Lyndcn-Bcll & Sandage 
(1962) (Sandage 1990, Majewski 1993, CMG97). The second infall episode, delayed to 
the first, forms the thin disk component, with a time-scale much longer than that of the 
thick disk. The infall rate /(r, t) is normalized to the local disk density at present time, 
i.e., f* g f(r, t) dt = £t t(r, t g ). Assuming the total masses of different rings in the disk at 
the present time are all exponentially decreased with the increasement of Galactocentric 
distance for thin and thick disk with the same scale- length r , the form of A{r) and B(r) 
in the pre-thin model can be written respectively as 

Air) = = — \ — e r a (6) 



l?» e v " dt 

J V27T0Y 

B(r) = < 



if t <C t max 

Stot(r ,t g )~S t fe»cfc(r0,t g ) r ,~ I -^ a - if +<- + <- + W 
I -It-Tj) 2 /2<r' i e L max _i L _i L g 

J tmax V 27rcr d 

respectively, where t g is the age of the Galactic disk, t max represents the epoch of time 
at which the formation of the thin disk begins, S tot (r ,t 9 ) is the present total surface 
density in the solar neighborhood, and ^thick{rQ,t g ) is the local surface density of the 
thick disk at the present time, which is one of the free parameters in our model. We 
adopted r = 2.7 kpc (Robin ct al 1996, Kent 1992), and E tot (r ,i g ) = 55.0 M Q p C - 2 
(Rana 1991, Sackett 1997) for Galactic disk. There are four free parameters in the pre- 
thin model, T t ,T d ,t max and T, thick (r Q ,tg). 

Contrary to the pre-thin model, the post-thin model assumes the first infall episode 
forms thin disk. Then, the thick disk forms delayedly as a result of some actions on or 
by the thin disk. Using the same method described above, the form of A(r) and B(r) in 
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the post-thin model can be written as 



if t < t 



max 



Mr) = { v th ^ k (ro,t g ) , (8) 

-(t-x d )2/2^ c 11 L max ^ (- ^ f-g, 



ft: 



■ dt 



B(r) = W?Q,t 3 )-^ ick(r Q ,t g ) iz^e 
Jo ^<j t m 

respectively, where each parameter has the same notation as Eq. (6) and (7) except t max , 
which represents the epoch of time at which the thick disk begins to form. There are also 
four free parameters in the pre-thin model, T t7 Td,t max and T, t hick(rQ,t g ). 

3.3. SFR and IMF 

In the majority of chemical evolution models, star formation rate (SFR) is assumed 
to depend on some power of the gas surface density ( Prantzos & Aubert 1995, Tosi 
1996, CMG97). Based on gravitational instability, Wang & Silk (1994) developed a self- 
consistent model to derive the global star formation rate as a function of radius in galactic 
disks. The resulted star formation rate not only depends on the gas surface density, but 
also is proportional to the epicycle frequency n. Since k oc r _1 , we adopt a similar star 
formation rate as that of Prantzos & Aubert (1995), which can be expressed as (in units 
of MQpc^Gyr- 1 ) : 

<P(r,t)=v^ as (r,t)/r (10) 

where S gas (r, t) and r are in units of M Q pc~ 2 and kpc, respectively. The power law index 
n = 1.4 is adopted (Prantzos & Silk 1998), which is in some degree similar to that of 
Kcnnicutt (1998). The value of v is derived from the condition of reproducing the present 
observed gas surface density in the solar neighborhood. We adopted S gas (r Q ,t g ) = 10.0 
M Q pc~ 2 (Scoville & Sanders 1987, Prantzos & Aubert 1995, Sackett 1997). 

The adopted stellar initial mass function (IMF) is taken from Kroupa et al. (1993), 
in which the IMF is described by a three-slope power law, <j)(m) cx m~( 1+x \ In the 
high-mass region, the IMF has a relatively steep slope of x = 1.7, while it flattens in 
the low-mass range (x = 1.2 for 0.5M Q < m < 1.0M Q and x=0.3 for m < O.5M ). The 
adopted IMF is normalized to /q" ' rruj){m) dm = 1. 

4. Results and discussions 

4.1. X 2 -test 

The observed G-dwarf metallicity distribution is treated as the tightest observational 
constraints on the chemical evolution models of the Galactic disk(CMG97). The metal- 
licity of local disk stars in the sample of Hou et al (1998) extends to [Fe/H] = -1.5, which 
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is consistent with the observations of the thick disk (Majewski 1993, Nissen & Schuster 
1997). However, the metal- weak tail of the thick disk is excluded from the distribution 
of RM96 by using a chemical criterion. Since the main aim of this work is to predict the 
general properties of the thick disk, we select the distribution of Hou et al (1998) as the 
observational constraint. 

The x 2 01 the goodness-of-fit to the observed data is calculated as follows 

n 

xw = 2^( )' ( n ) 

»=i ai 

where y m i and y i are respectively the model-generated and observed data in ith data 
bin, ai is the error of the observed data in ith bin, and n is the total number of bins. 
For the metallicity distribution, Hou et al (1998) have already divided the observed data 
into 19 bins. The error Oi is given by 

,AiV \ JKNl f ANj 1 ~ IAN, 1 AiV, 1 

"hr* = v ( ^^ } + (_ r ^ } = V + hr^ ' (12) 

where ANi/N is the relative number of G-dwarfs in the ith bin and N is the total number 
of G-dwarfs in the corrected sample. 

For comparison, we consider four different models respectively, which are closed-box, 
one-component, pre-thin and post-thin model. The differences among these models are 
the treatments of the infall rate. We perform model calculations for a broad range of free 
parameter combinations, if there are any, to present the best-fit to the new metallicity 
distribution. Using x 2 -test, the results of best-fit to the G-dwarf metallicity distribu- 
tion for these individual models are shown in table 2, which are detailedly discussed at 
following subsections. The characteristics of the models are shown in the first column. 
The resulted values of free parameters, the \ 2 an d the model confidential levels for the 
best-fit are shown in second, third, fourth column respectively. Figure 1 presents the best- 
fits to the new G-dwarf metallicity distribution for closed-box model (dotted line), the 
one-component model (dot-dashed line), pre-thin model (long dashed line) and post-thin 
model (full line). 

4-2. Closed-box model 

The closed-box model considers the whole disk as an isolated system. Thus, A(r) = 
B(r) = 0. No free parameter exists if the closed-box model is considered. Figure 2 shows 
that the closed-box model predicts larger number of metal-poor stars than that of the 
observations. It is often called the G-dwarf problem. The fact that closed-box model does 
not work is also shown in table 2, where the confidential level is too low to accept the 
model. Several possible explanations to the G-dwarf problem have been proposed (see 
brief reviews in Francois et al 1990, Malinie et al 1993). Generally, a good agreement 
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Table 2. The results of best-fit to the new G-dwarf metallicity distribution for four 
different models 



characteristics of 


free parameters 


x 2 


model confidential level 


the model 


(best- fit) 


(best-fit) 


(best-fit) 


closed-box model: 








A(r) = 








B(r) = 


no 


131.7 


< 0.1% 


one-component model: 








A(r) = 








B(r) # 


r d = 3.8Gyr 


20.3 


30% 


pre-thin model: 








Ayr) 7^ 


Tt=l.O Gyr 






B(r) / (=0, if t < tmax) 


r d =4.0Gyr 


14.6 


70% 


tmax 7^ 


£ tWcfe (r Q ,t 9 ) = 14.0M Q pc- 2 
t max = l.QGyT 






post-thin model: 








A(r) / (=0, if t < tmax) 


r t =1.0 Gyr 






B(r) + 


r d =4.0Gyr 


14.4 


70% 




£ thlcfc (r ,t 9 ) = lO.OM Q pc- 2 
t m ax = 1.0Gyr 







between model predictions and observations is obtained by models that assume the disk 
formed by infall of piromordial gas (Pilyugin & Edmunds 1996b, CMG97) 

4-3. One-component model 

The one-component model treats the thick and thin disk as one disk, which corresponds 
the infall rate of only one Gaussian form, i.e., A(r) = 0, B(r) ^ 0. There is only one free 
parameter t^, i.e., the infall timescale of the whole disk in one-component model after 
normalization mentioned above. Figure 1 shows that , for the one-component model, 
the goodness-of-fit between best-fit model predictions and observations seems to be ac- 
ceptable, especially for the shape of the single peak of the metallicity distribution for 
G dwarfs. But, the results of x 2 -test gives that the confidential level of one-component 
model is some 30%, yet it can reach 70% for the two-component model (Table 2). This 
means one-component model is not the best, which is consistent with the observational 
result that the thick disk is kincmatically and chemically different from the thin disk (see 
brief reviews in Majewski 1993). Therefore, it is necessary to treat the thick and thin 
disk differently, if one want to investigate detailedly the Galactic chemical evolution. 
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Fig. 1. The best- fits to the new G-dwarf metallicity distribution for close-box model 
(dotted line), the one-component model (dot-dashed line ), pre-thin model (long dashed 
line) and post-thin model (full line). The observed data are taken from Hou et al (1998). 

4-4- Two-component model 

For the two-component model, both the pre-thin and post-thin model are considered. 
The differences between these two models are represented by different treatments of the 
infall rate, which are detailedly described in section 2.2. There are four free parameters in 
both the pre-thin and post-thin model, Tt, r^, t ma x, ^thick{fQ, t g ). It should be emphasized 
that the t ma x has different meaning in different model. Figure 1 show that, for these two 
models, the best-fit model predictions are in good agreement with the observations, which 
are confirmed by the results of quantitative tests (Table 2). From Figure 1 and Table 2, 
it is difficult to distinguish which model is better. 

One important parameter predicted by the best-fit of our two-component model is 
the local surface density of thick disk at the present time. On the other hand, the density 
ratio of thick to thin disk are usually deduced from studies of star counts. In Table 3, 
we present the comparison between our model predictions and the data compiled from 
literatures. Since the density normalization and scale-height of the thick disk are anti- 
correlated when fitted simultaneously (Reid & Majeweski 1993, Majeweski 1993, Robin et 
al. 1996), the scale heights of the thick disk are also shown in 2th column. The previous 
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Table 3. Local density ratio of thick to thin disk from literatures 



Authors 


scale height of thick disk 


space density 


^thick(r Q,tg) 




(PC) 


normalization 


M e pc~ 2 


Gilmore & Reid (1983) 


1450 


0.020 


4.85 


Rose (1985) 


500 - 1000 


0.100 


11.0 


Friel (1987) 


1000 


0.050 


7.85 


Sandage (1987) 


940 


0.110 


14.1 


Norris (1987) 


1100 


0.034 


6.09 


Kuijken & Gilmore (1989) 


1000 


0.040 


6.47 


Reid & Majewski (1993) 


-i a r\r\ i f r\r\ 

1400 ~ 1600 


0.020 ~ 0.025 


5.67 


Robin et al. (1996) 


760 


0.056 ± 0.010 


6.83 


Ojha et al. (1996) 


760 


0.074 ± 0.020 


8.69 


Buser et al. (1998) 


1150 


0.054 ± 0.015 


9.43 


This paper (1999): 








pre-thin model 


1000 


0.113 


15.0 


post-thin model 


1000 


0.067 


10.0 



results for the space density ratio of thick to thin disk (the parameter D) are shown 
in 3th column. Taking the scale-height of thick and thin disk as h t — lOOOpc and hd = 
300pc respectively as usual (Majewski 1993), the local space density ratio of thick to thin 
disk at the present time predicted by our model can be obtained based on the following 
equation: 

p ^ ^thick(rQ,t g ) ha 

^tot(rQ,t g ) — Y,thick(rQ,t g ) h t 

where £t/u C fc (?"© i i g )=10.0M Q pc~ 2 and 15.0 M & pc~ 2 for the best-fittings of the pre-thin 

and post-thin model, respectively. Moreover, equation (13) can also be used for deducing 

the Ti t hick{rQitg) with published D and ht in literatures. The results are shown in 4th 

column of Table 3 with ft,£(=300pc and E tot (r0, t g ) =55. 0Mqpc~ 2 , respectively. 

From Table 3, it can be seen that the previous density normalizations of thick disk 

span a wide range (from 0.02 to 0.11). This is partly due to the difficulty in distinguishing 

the kinematical, chemical and spatial characteristics of thick disk with halo and thin 

disk. Moreover, comparing published star counts and color distribution in neighboring 

fields, star density discrepancies are sometimes larger than photometric random errors 

as established by authors (Ojha et al. 1996, Robin et al. 1996). But the results of more 

recent surveys are in good agreement within the error range (Ojha et al. 1996, Robin 

et al. 1996, Buser et al 1998). Table 3 shows that the value of £t/»icfc (?"© ; t g ) predicted 

by the post-thin model are consistent with most of the previous results from literatures, 

while ^thick{ r Qj tg) predicted by the pre-thin model is larger than the previous data 
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Fig. 2. The \ 2 as a function of the 'SthickifQjtg) within reasonable range for both the 
pre-thin model (long dashed curve) and post-thin model (full curve) in which other free 
parameters are fixed as r t = 1.0Gyr, t TOaa ,=1.0Gyr and T c ;=4.0Gyr. Three short horizontal 
lines indicate the model confidential level of 30% (dotted line), 50% (long dashed line) 
and 70% (full line), respectively. The points indicate the different results of x 2 if the 
value of ^thick(re,t g ) from ten literatures (Table 3) are adopted for both the pre-thin 
model (open squares) and post-thin model (open circles), respectively. 

from studies of star counts. To illustrate this quantitatively, Figure 2 shows the x 2 as a 
function of the Ti t hick{fQ,tg) within reasonable range for both the pre-thin model (long 
dashed curve) and post-thin model (full curve) in which other free parameters arc fixed as 
Tt=1.0Gyr, t max =1.0Gyr and Td—^-OGyr. It is shown that the value of x 2 is very sensitive 
to Ti t hick{fQn tg)- This suggests that the thick disk has great influence on the Galactic 
chemical evolution. In figure 2, three short horizontal lines indicate the model confidential 
level of 30% (dotted line), 50% (long dashed line) and 70% (full line), respectively. The 
points in Figure 2 indicate the different results of \ 2 if t ne value of ^thicki^Q^tg) from 
ten literatures (Table 3 ) are adopted for both the pre-thin model (open squares) and 
post-thin model (open circles), respectively. 

Figure 2 shows that, for the post-thin model, five points have model confidential level 
larger than 70%, while only two points have confidential level larger than 70% for the pre- 
thin model. This suggests that the post-thin model be better than pre-thin model. Other 
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Fig. 3. Temporal evolution of [Fe/H] predicted by the best-fit model (full line) for the 
solar vicinity. The observed data (asterisks) are from Edvardsson et al (1993). 

evidences tended to favour the post-thin scenario for the formation of thick disk comes 
from the following facts. First, the thick disk is kinematically distinct from the thin disk 
and it shows no kinematic gradients (Ojha et al 1994 a,b). Second, Gilmore et al (1995) 
studied metallicity distribution of thick disk stars up to about 3kpc from the Galactic 
plane. They found that thick disk stars show no vertical abundance gradient. This argues 
against dissipational setting as the formation process of the thick disk (Freeman 1996). 

4-5. Discussions 

Based on the above discussions, we treat the post-thin model with rt=1.0 Gyr, 
imax = l-0Gyr, 7^=4.0 Gyr and ^thick{ro,t g )= 10.0Mqpc~ 2 as the best-fit model. 

Figure 3 presents the comparison between our best-fit model predictions for the AMR 
and the observations. The full line is our model predictions and the points are observa- 
tional data taken from Edvardsson et al (1993). Figure 3 shows that, at the beginning 
of the formation of the thick disk [t = t max = l.OGyr), the iron abundance of ISM de- 
creases a little due to the increasing infall rate of the primordial gas. After that phase, the 
model predicts that the metallicity increases smoothly with time. The overall tendency 
for this relation is consistent with the mean observations, but the present model can not 
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Fig. 4. Predicted behavior of [O/Fe] vs. [Fe/H] for the best fit model (full line). The 
observed data are from Edvardsson ct al (1993) (asterisks) and Barbuy (1988) (full 
squares) . 

reproduce the large observed scatters. Nordstrom et al. (1997) discussed in detail the 
main hypotheses for the origin of this scatter, such as star formation in an inhomoge- 
neous gaseous medium, orbital diffusion in homogeneous galaxy and mergers or accretion 
events. However, a physical mechanism that reproduces the observed scatter in the AMR 
without violating other observational constraints has not yet been identified. 

Figure 4 compares the predicted behaviors of [O/Fe] vs. [Fe/H] for the best-fit model 
(full line) with the observations. The observed data are taken from Edvardsson et al 
(1993) (asterisks) and Barbuy (1988) (full squares). Our model predicts that there is 
a small loop at [Fe/H]= -1.6. The similar behavior is predicted in the best-fit model 
of CMG97. Figure 4 shows that our model prediction is in good agreement with the 
observations. This suggests that the relative stellar yield for oxygen and iron we adopted 
here be reasonable. 

Contrary to the case of the solar neighbourhood, the available observations for the 
Milky Way disk offer information maily about its current status, not its past history. 
Therefore, there is much more freedom in constructing a model. Up to now, chemical 
evolution models of the Galactic disk consider the disk as a system of independent rings. 
This oversimplification generally ignores the possibility of radial inflows produced in 
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Fig. 5. Radial distribution of the oxygen abundance at the present time predicted by the 
best-fit model (full line). The observed data of HII regions are from Shaver et al (1983) 
(asterisks), Fich & Silkey (1991) (open circles) and Vilchez & Esteban (1996) (open 
squares) . The observations of early B-type main-squence stars are taken from Smartt et 
al (1997) (full squares) and Gummersbach et al (1998) (crosses). 

gaseous disks, e.g. by viscosity or by the infall of gas with a specific angular momentum 
different from that of the underlying disk (Prantzos & Silk 1998). Fortunately, a radial 
variation of the infall time-scale may play a similar role. In our best-fit model, the infall 
timescale of the thin disk is assumed to be radially dependent, taking an lower values 
in the inner disk (t<2=2 Gyr at r—2 kpc) and larger ones in the outer disk (rd—A Gyr at 
r=8.5 kpc). 

In Figure 5, we compare the predicted radial distribution of oxygen abundances by our 
best-fit model (full line) with the observations. The observed data of HII regions are from 
Shaver et al (1983) (asterisks), Fich & Silkey (1991) (open circles) and Vilchez & Esteban 
(1996) (open squares). The observations of early B-type main-squence stars are taken 
from Smartt et al (1997) (full squares) and Gummersbach et al (1998) (crosses). It shows 
that model predictions of oxygen abundances are larger than most of the observations. 
Moreover, our model predicts that the abundance gradient in inner region (-0.01 dex/kpc 
for r < 8.5 kpc) is steeper than that in outer region (-0.07 dex/kpc for r > 8.5 kpc). 
This seems to be contradictory to the observations of HII region, which suggests a flatter 
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oxygen abundance gradients in outer region. This disagreement was also obtained by 
most of chemical evolution models for the Galactic disk (see brief reviews in Tosi 1996). 
This is probably due to the simplicities of the present models. Samland et al. (1997) 
developed a chemi-dynamical evolution model for the Galactic disk and presented better 
fit for the observed variations of oxygen abundance across the whole disk. Therefore, 
this disagreement probably could be solved if one considers the influence of dynamical 
evolution on the chemical evolution and the effect of gas heating on star formation rate. 

Figure 6 presents the comparison between the predictions of our best-fit model for 
the radial profile of the present SFR (full line) and the observations. Observed data 
are normalized to SFR in the solar neighborhood. Data are based on several tracers of 
star formation: Lyman continuum photons from HII regions (full squares, from Gusten 
& Merger 1983); pulsars (open circles: from Lync et al 1985); and supernova remnants 
(crosses: from Guibert et al 1978). The good agreement between our model prediction 
and the observations indicates that our model predicts reasonable star formation history 
for the Galactic disk. 

Finally, in Figure 7, we present the comparison for the radial distribution of the 
present gas surface density between our best-fit model predictions (full line) and the 
observations (dashed lines). Two dashed lines are reproduced from Prantzos & Aubert 
(1995). The lower dashed line is the sum of atomic and molecular hydrogen given in 
Dame (1993), corrected for the contribution of 30% helium. The upper one is obtained 
by adopting the gas surface density in the solar neighbourhood as 16 M Q pc 2 and scaling 
the curve of Dame (1993) accordingly (Prantzos & Aubert 1995). Given the uncertainties 
in the observational data, the model is also in good agreement with the observed profile. 

5. Conclusions 

In this work, we introduce a two-component models for the chemical evolution of the 
Galactic disk, which assumes that the formation of the thick and thin disks occur in two 
main accretion episodes. The infall rate is assumed to be Gaussian. Both the pre-thin and 
post-thin scenarios for the formation of the Galactic disk are considered. The local surface 
density of the thick disk at the present time is chosen to be one of the free parameters. 
Following Prantzos & Silk (1998), we also assume that the SFR is not only proportional 
to n power of gas surface density, but directly correlates with the Galactocentric dis- 
tance. Comparing model predictions with the new metallicity distribution in the solar 
neighbourhood (Hou et al 1998), we use the x 2 ~test to derive best-fittings and com- 
pare the reasonableness of four different models, which are closed-box, one-component, 
pre-thin and post-thin models. Moreover, comparisons between the predictions of our 
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Fig. 6. Profile of present star formation rate (SFR) in the galactic disk. Full line is the 
predictions of the best-fit model. Observed data are normalized to SFR in the solar 
neighborhood. Data are based on several tracers of star formation: Lyman continuum 
photons from HII regions (full squares, from Gusten & Merger 1983); pulsars (open 
circles: from Lyne et al 1985); and supernova remnants (crosses: from Guibert et al 
1978). 

best-fit model and the main observational constraints are presented. The main results 
can be summarized as follows. 

1. Our results suggests that the post-thin model for the formation of the Galactic disk 
should be preferred. This is consistent with the observational evidences that the thick 
disk is chemically and kinematically distinct from the thin disk and it shows no 
vertical abundance gradient. 

2. The goodncss-of-fit for model predictions about metallicity distribution to the obser- 
vations, x 2 , is very sensitive to the local surface density of the thick disk at the present 
time. This suggests that it is necessary to treat the thick and thin disks differently if 
one want to investigate detailedly the Galactic chemical evolution. 

3. The post-thin model predicts ^thick{ r Qi t g ) = IO.OMqpc^ 2 . The resulted space density 
ratio of thick to thin disk is consistent with the previous data from recent studies of 
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Fig. 7. Radial distribution of the present gas surface density predicted by the best-fit 
model (full line). The lower dashed line is the sum of atomic and molecular hydrogen 
given in Dame (1993), corrected for the contribution of 30% helium. The upper one is 
obtained by adopting the gas surface density in the solar neighbourhood as 16 Mqpc 2 
and scaling the curve of Dame (1993) accordingly (Prantzos & Aubert 1995). 

star counts. However, the pre-thin model predicts a larger value of the local thick 
disk density. 

4. The predictions of our best-fit model are in good agreement not only with the observed 
data in the solar neighborhood, but also with the main observational features of the 
Galactic disk. However, contrary to the observations in HII regions, our model predicts 
the oxygen abundance gradient in outer region is steeper than that in inner region. 
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